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Abstract 

A method  of  introduction  of  the  effective  material  tensors  of  bianisotropic  multilayered  peri- 
odic structures  based  on  the  approximate  calculation  of  the  characteristic  matrix  of  the  unit 
cell  of  the  system  with  the  help  of  Campbell-Hausdorff  series  [1]  is  developed.  Obtained 
effective  tensors  are  valid  for  the  use  in  a wide  wave  band.  This  paper  is  primarily  concerned 
with  a comparison  of  the  accuracy  of  different  approximations,  namely,  a long  wavelength 
approximation  and  approximations,  using  three  and  five  terms  of  Campbell-Hausdorff  series. 

1.  Introduction 

Attention  now  frequently  focuses  on  the  study  of  the  effective  properties  of  plane  stratified 
periodic  media.  In  doing  so  most  often  a long  wavelength  approximation  is  used.  For  example, 
in  [2]  the  coordinate-free  formulae  for  the  effective  tensors  of  permittivity  e,  permeability  \i  and 
pseudotensors  of  gyrotropy  a,  /3  of  the  plane  stratified  periodic  bianisotropic  systems  of  the 
most  general  type  were  obtained.  But  it  is  well  known  [3],  [4]  that  some  systems  composed 
from  nongyrotropic  layers  can  possess  gyrotropic  properties  due  to  their  specific  structure  (i.e. 
possess  form  gyrotropy).  Phenomena  of  this  kind  can’t  be  explained  in  the  framework  of  a long 
wavelength  approximation.  To  treat  such  problems  it  is  necessary  to  extend  the  wave  band,  in 
which  the  effective  tensors  can  be  used.  One  way  of  doing  this  is  to  employ  an  approach  [5], 
based  on  the  approximate  calculation  of  the  characteristic  matrix  of  the  unit  cell  with  the  help  of 
Campbell-Hausdorff  series.  Recently  it  was  shown  that  by  using  this  approach  the  effects  of  form 
gyrotropy  (bianisotropy)  can  be  really  explained  in  terms  of  the  theory  of  effective  parameters 
(see  [6],  for  example).  In  [6],  [7]  our  consideration  was  limited  by  the  three  terms  of  this  series. 

Here  we  obtain  more  accurate  formulas  using  the  forth  and  fifth  terms  of  Campbell-Hausdorff 
series. 

We  consider  systems  formed  by  a periodic  set  of  plane  bianisotropic  layers  with  different 
thicknesses  ln  (n  = 1, 2, ...  N,  where  N is  the  number  of  the  layers  forming  the  unit  cell).  The 
layers  are  characterized  by  the  constitutive  relations 

Dn  = "t"  flnHn,  Bji  = "h  (1) 

where  en  ,/zn  and  an,  /3n  are  the  dielectric  permittivity,  the  magnetic  permeability  tensors  and 
the  pseudotensors  of  gyrotropy,  respectively. 

2.  Effective  Material  Tensor  Parameters 

The  characteristic  matrix  P = exp  ( ikolM ) of  a layer  with  thickness  l relates  the  six-vectors 
(E,H)t  at  the  layer  boundaries  (ko  = lo/c).  M is  some  matrix  depending  on  the  parame- 
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ters  of  the  layers,  angle  and  plane  of  incidence.  Matrix  M is  given  in  an  explicit  form  in  [7]. 
Henceforward  we  shall  follow  notation  used  in  [7]. 

Let  us  consider  the  system  formed  by  two  alternate  layers  with  different  thicknesses  ln  and 
different  sets  of  tensor  constants  en,  /in,  an,  /3n  (n—  1,2).  In  this  instance  the  characteristic 
matrix  of  the  unit  cell  has  the  form 


P — exp  ( ikoLM ) = exp  (zfco/2M2)  exp  ( ) , 


(2) 


where  L = l\  + l2  is  the  system  period  and  M is  the  matrix  to  be  found.  This  matrix  can 
be  expressed  in  terms  of  the  layers  parameters  with  the  help  of  Campbell-Hausdorff  series  as 
follows 

M = /iMi  + /2M2  + ifi  ^ [M2,  Mi]  - P*2,  Mi] , M\]  + f2  [[. MUM2 ] , M2]}  + • • • (3) 

where  [M2,Mi]  = M2M\  — MiM2.  If  i dn/X  <C  1,  then  the  series  (3)  quickly  converges  and  one 
can  drop  the  remainder  of  it  after  the  fc-th  term.  The  first  and  the  second  terms  correspond  to 
the  long  wavelength  limit.  The  possibilities  of  introduction  of  the  effective  constitutive  tensors 
having  regard  the  third  term  are  discussed  in  [7].  In  what  follows  we  shall  break  the  series  after 
the  fifth  term. 

In  perfect  analogy  with  [7]  we  shall  use  in  (2)-(3)  the  matrices  M/,Mi/,M2/,  describing 
the  transformation  of  the  tangential  components  of  the  field  vectors,  instead  of  M,M\,M2, 
and  we  assume  q to  be  the  left  and  the  right  eigenvector  of  each  of  the  tensors  en,  nn,  «n,  (3n 
(£nq  = q£n  = £ngq,-*-)- 

To  compare  the  accuracy  of  different  approximations,  namely,  long  wavelength  approxima- 
tion and  approximations,  using  three  and  five  terms  of  Campbell-Hausdorff  series,  let  us  consider 
the  most  simple  case  of  normal  incidence  onto  the  nongyrotropic  nonabsorbing  layers  of  the  same 
thickness  l\  = l2  = l.  As  this  takes  place  (3)  reduces  to 

Mj  = - (Mi/  + M2I)  + i— {[M2j,  Mu]  + [[M2/,  Mu] , Mu  - M2/]}  -I (4) 

Generally  speaking,  Campbell-Hausdorff  series 

Z = X + Y + ±[X,Y]  + ±[X,[X,Y)]  + -^[Y,[Y,X}}  + - (5) 


converges  provided  ||X||  < ||Y||  < 

of  X.  In  our  case  series  converges  if 


where  ||A||  = 


trace 


is  Euclidean  valuation 


7rmax(li,l2 ) ln2 

X ~ 4 max(\\Mu\\,  \\M2T\\) 


(6) 


For  example,  in  the  case  of  nonmagnetic  crystals  with  ||e/||  ~ 4.5  (eigenvalues  of  £j  ~ 3.0) 
^ < 0.039. 

After  simple  calculation  we  find  the  effective  material  tensors  of  the  system  at  hand 


+ 7 xl 

<*i  = Pj  = *2X^0’  Ra  = (£l/qX - ^£/qxMi/)  > 
£J  = \ (£1^  + £2i)  + e\  VI  = \ (/^u  + m)  + v', 
e'  = i ^ R'e , K = -i?Qqx  Ae/  - Ae/q><^, 


(7) 

(8) 
(9) 
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/*'  = R'n,  R'n  = Rpcf&Vr  + A^/qxi?a,  (10) 

where  A e/  = £21  - £1 7,  A m = /z2/  - 

It  is  worth  noting  that  third  term  of  Campbell-HausdorfF series  give  rise  to  the  effective  ’’ten- 
sors of  gyrotropy”  aj,  Pi  while  forth  and  fifth  terms  contribute  only  in  the  effective  permittivity 
and  permeability  tensors  (terms  e',  yJ  in  (8)). 

^From  (7), (9),  (10)  it  follows  that 


n*j<2(iiei,ii  iiawIi  + iia£/||  iif,, id, 

(11) 

M<2(|M|  ||A/x,||  + ||Ae,||  Ilf „ II), 

(12) 

Ml  < 2 IIAejH  (pyi  + ||Bf  ||) , 

(13) 

KII<2||Af,||(||ii0||  + ||iJs||). 

(14) 

^From  above  estimations  and  formulae  (7)- (10)  it  is  clear  that  in  the  convergence  range  the 
contribution  of  the  forth  and  fifth  terms  of  Campbell-Hausdorff  series  in  effective  tensors  is  far 
less  than  the  contribution  of  third  term,  and  usually  may  be  thought  of  as  negligibly  small  if 
compared  with  the  contribution  of  the  first  and  second  terms.  But  if  in  the  long  wavelength 
approximation  the  system  under  hand  has  transversely  isotropic  permittivity  and  permeability 
tensors,  then  contribution  of  the  forth  and  fifth  terms  can  be  noticeable.  The  latter  case  is  of 
great  concern  in  investigation  of  form  bianisotropy  (gyrotropy),  because  the  manifestations  of 
gyrotropy  are  often  suppressed  due  to  the  permeability(permittivity)  being  tensors,  especially 
in  the  optical  wave  region  [9]. 

3.  Conclusion 

In  this  paper  we  discussed  primarily  the  comparative  contribution  of  different  terms  of 
Campbell-Hausdorff  series  in  the  effective  tensors  of  ’’permittivity,  permeability”  and  ’’gy- 
rotropy” , analyzing  the  simple  case  of  normal  incidence  onto  the  system  of  nongyrotropic  layers. 
Obtained  estimations  holds  true  for  gyrotropic  layers  too,  at  least  for  the  short  enough  wave- 
lengths. At  normal  incidence  generalization  of  formulae  (7)-(10)  to  the  case  of  bianisotropic 
layers  is  not  a particular  problem. 

Of  course,  derived  effective  tensors  are  not  true  constitutive  tensors  because  generally  they 
are  valid  only  in  the  immediate  vicinity  of  normal  incidence.  It  is  possible  to  introduce  true 
constitutive  tensors  (i.e.  not  depending  on  the  angle  and  plane  of  incidence)  even  having  regard 
to  five  terms  in  Campbell-Hausdorff  series,  but  only  for  some  systems  with  the  specific  relation 
between  the  layers  parameters.  Therefore,  proposed  method  is  more  convenient  for  analysis 
of  the  effective  properties  of  systems  at  normal  incidence,  when  there  is  no  restrictions  on  the 
parameters  of  the  layers.  Extensive  computations  made  with  the  use  of  exact  and  approximate 
formulas  are  in  good  accordance. 
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